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Abstract 

Partial agonists exhibit a submaximal capacity to enhance the coupling of one receptor to an intracellular binding partner. 
Although a multitude of studies have reported different ligand-specific confornnations for a given receptor, little is known 
about the mechanism by which different receptor conformations are connected to the capacity to activate the coupling to 
G-proteins. We have now performed molecular-dynamics simulations employing our recently described active-state 
homology model of the dopamine D2 receptor-Ga, protein-complex coupled to the partial agonists aripiprazole and 
FAUC350, in order to understand the structural determinants of partial agonism better. We have compared our findings 
with our model of the D2R-Gai-complex in the presence of the full agonist dopamine. The two partial agonists are capable 
of inducing different conformations of important structural motifs, including the extracellular loop regions, the binding 
pocket and, in particular, intracellular G-protein-binding domains. As G-protein-coupling to certain intracellular epitopes of 
the receptor is considered the key step of allosterically triggered nucleotide-exchange, it is tempting to assume that 
impaired coupling between the receptor and the G-protein caused by distinct ligand-specific conformations is a major 
determinant of partial agonist efficacy. 
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Introduction 

G protein-coupled receptors (GPCRs) constitute an important 
class of membrane-bound glycoproteins that participate in the 
regulation of various physiological processes including heartbeat, 
breathing and our senses of vision, smell and taste [1]. By 
transmitting extracellular stimuli to the inside of the cell, GPCRs 
serve as linker molecules that connect ligand binding to the 
coupling of intracellular binding partners including G proteins or 
P-arrestins [2]. In its canonical signaling pathway, the activated, 
neurotransmitter/hormone-occupied receptor couples to G pro- 
teins, thereby inducing conformational changes that give rise to 
nucleotide-exchange and culminate in various functional responses 
[3-6]. The efficacy of a given ligand refers to its capacity to 
enhance the coupling of one receptor to a particular intracellular 
effector protein, thereby inducing quantifiable cellular responses 
[7,8]. Depending on the extent of their functional response, 
ligands can be classified as neutral antagonists or inverse agonists, 
neither of which stimulates receptor activation, fuU agonists, which 
cause a cellular response that strongly resembles that of the 
endogenous ligand, and partial agonists, which exliibit submax- 
imal effects even at saturating concentrations. 

Partial agonism at dopamine D2 receptors (D2R) has been 
suggested to exert beneficial effects on schizophrenia, a chronic 
mental illness characterized by hypo- and hyperfunctions in 
monoamine neurotransmitter systems, including the mesolimbic 



and mesocortical dopaminergic pathways [9]. The dopamine 
receptor partial agonists aripiprazole and the drug candidate 
cariprazine represent promising options for the treatment of 
schizophrenia [10-12] because of their stabilizing effect on 
monoamine pathways, especially the dopaminergic pathways, 
and their atypical antipsychotic effect. 

Understanding the molecular basis of partial agonism requires 
detailed insight into the impact of ligands with different efficacies 
on the conformation of a given receptor or receptor-effector 
complex. Indeed, earlier studies that explored the origin of partial 
agonism revealed a ligand-dependent modulation of micro- 
switches important for receptor activation [13,14] and specific 
ligand-specific conformations within receptor epitopes that include 
the orthosteric binding pocket, the extracellular loop (EL) region 
and areas of receptor-G protein coupling [14-20]. However, a 
structural study based on an active-state receptor coupled to an 
intracellular binding partner able to link ligand-specific receptor 
conformations to the capacity to activate a given effector has not 
yet been reported. 

Significant progress in the field of GPCR-crystaUography led 
the way to the elucidation of the crystal structure of the ternary P2 
adrenergic receptor G^ protein complex [21], which provides an 
unprecedented framework for the investigation of ligand-induced 
receptor conformations and of receptor G protein interactions. 
Using this experimental structure as a template, we were able to 
generate active-state homology models of the dopamine D2 
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receptor (D2R) in complex with Ga; and the endogenous agonist 
dopamine. We used these models to explore the structural 
determinants of selective receptor-G protein coupling on the 
amino-acid level using microsecond molecular-dynamics (MD) 
simulations [22]. 

We now report long-term all- atom MD simulation.s of our D2R- 
G(Xi homology model designed to help understand the molecular 
basis of partial agonist activity better. The D2R-Go(i model in a 
hydrated bilayer environment was coupled to the partial agonist 
aripiprazole and a structurally related 1,4-disubstituted aromatic 
piperazine (1,4-DAP) [23,24], FAUC350, which has been shown 
to induce a strongly impaired modulation of cAMP accumulation 
while activating ERKl/2 phosphorylation with moderate efficacy 
[25]. In this work, we have investigated the impact of aripiprazole 
and FAUC350 on the conformation of the ternary D2R-GQ(i- 
complex, focusing on the shape of the extracellular loop region, 
the binding pocket and receptor-G protein-binding epitopes and 
compared our findings with our homology model of the D2R-Gai- 
complex in presence of the full agonist dopamine. As in our 
previous work [22,26,27], we have performed single long 
simulations, rather than multiple shorter ones. This strategy 
allows us to avoid missing conformational changes that occur with 
a characteristic induction period, which either may be inherent to 
the natural system or be caused by a parameterized kinetic or 
thermodynamic structural bias of the force field towards "native" 
(i.e. X-ray or homology) protein structures. 

Results/Discussion 

Molecular-dynamics Simulations of Ternary D2R-Gai- 
complexes Coupled to Dopamine, Aripiprazole and 
FAUC350 

MD simulations with the ternarv- dopamine D2R-Gai complex 
revealed a high degree of conformational fl('xibility for the x\- 
angle of His393 (atoms: C-Ca-Cfi-Cy) [22], a crucial residue 
for D2R activation [25,28]. Despite this flexibility, a dihedral 
angle of around —60°, which was found almost continuously in the 
simulation denoted D2'^''R (Table SI), turned out to be connected 
to favorable ligand-receptor, receptor-receptor and receptor-G 
protein-interactions (Figure SI). This simulation was used as a 
reference model with which to compare the simulations of D2R- 
Ga; coupled to the partial agonists aripiprazole and FAUC350 
(Tables SI and S2). Starting from the existing membrane-inserted 
model, the dopamine ligand was removed and the partial agonists 
were docked. A more detailed description is given in the Materials 
and Methods section. The aripiprazole- and FAUC350-bound 
D2R-G0!, complexes were subsequently submitted to energy 
minimization, equilibration and MD simulations runs of 800 ns 
and 500 ns, respectively. The G protein-moieties were more 
mobile than the receptor-units in both simulation systems (Figures 
82 and S3). Nevertheless, the complexes showed no tendency for 
the G protein to dissociate from the receptor throughout the 
simulations (Figure S4). 

Aripiprazole and FAUC350 Open the Binding Pocket 
towards the Extracellular Surface 

The models show that dopamine and the phenylpiperazine 
moieties of aripiprazole and FAUC350 occupy the same 
orthosteric binding pocket, thereby interacting with residues of 
transmembrane helices (TMs) 3, 5, 6, 7 and EL2. The hnker 
moieties of the partial agonists show additional interactions with 
residues located at an extended binding pocket closer to the 
extracellular surface of the receptor (Figures S5 and S6). 
Depending on the Ugand, we observed differently shaped 



structures above the binding pocket (Figure 1), which appeared 
to close around the full agonist dopamine so that EL2 and the 
extracellular tail of TM7 (measured as the distance between 

residues Ilel83^'"^ and Tyr^ '^'^, Figure S7) are close to each other. 
In contrast, the partial agonists aripiprazole and, even more, 
FAUC350 opened up the binding pocket to the extracellular 
surface (Figure IB). This opening is connected with the formation 
of a second binding pose for the partial agonist FAUC350 (Figures 
S2 and S5). Our observations are consistent with previous studies, 
which suggested ligand-specific conformations for EL2 and a 
regulatory function for this loop in receptor acti\'ation [29,30]. 
The pronounced movement of the extracellular tail of TM7 
towards EL2 within the dopamine-complex can be attributed to a 
hgand-induced hydrogen bond between Ser^ '^^ and His^^^, which 
shifts the side chain of His®'^^ in the direction of TMS and clears a 
space for the inward movement of Tyr' "''^ (Figure 1). Unlike 
dopamine, aripiprazole and FAUC350 lead to a different dihedral 
angle (approximately 60° and, mainly, 180°) for His*' "^'^ (Figures lA 
and 2A), thereby increasing the distance between Tyr' ''* and EL2. 

Ligand-specific Interhelical Networks 

Residue His*' "^'^ was earlier found to play an important role in 
binding 1,4-DAPs, receptor activation and biased signaling 
[24,25,28]. The ligands dopamine, aripiprazole and FAUC350 
can induce different conformations within the binding pocket of 
D2R by specifically adjusting the dihedral angle of His*'"^'' 
(Figure 2A). Thus, we attribute a key role to His^"" in the 
ligand-dependent regulation of the binding pocket, which is 
consistent with experimental and computational reports. In 
addition to the ligand-specific conformationcd changes of His®'^^, 
we captured difierences in the interactions between the ligands 
and Ser^ *^, Ser^ *^ and Ser^ *^ in TM5 (Figure S5). These residues 
have been shown to be crucial for the binding of different ligands, 
including catecholamines, and for an effective receptor-G protein- 
couphng [31,32]. The full agonist dopamine formed stable 
hydrogen bonds to Ser'^ *^ and Ser^ *^ and stabilized a conforma- 
tion of Ser^ *^, which facilitated hydrogen bonding to His®'^^, 
while both aripiprazole and FAUC350 lacked hydrogen bonds to 
either serine residue and prevented the interaction between Ser'^ *^ 
and His^'^^ (Figures 1 and S5). Moreover, the frequency with 
which the side chain of Ser^ *^ pointed into the binding pocket was 
found to be reduced for aripiprazole and FAUC350 (Figure 2B). 
As receptor activation has been shown to be accompanied by such 
an inward movement of TM5-serines [33,34], the hindered 
engagement of these amino acids observed for aripiprazole and 
FAUC350 is likely to result in reduced a(:ti\'ation. In general, the 
ligands investigated exerted different effects on specific conforma- 
tions of amino-acid networks within the orthosteric binding pocket 
of D2R. Whereas the full agonist dopamine stabilized interactions 
between TMS and TM6 close to the binding pocket, aripiprazole 
and FAUC3S0 led to helical reorientations su[:h that TM6 moved 
away from TMS towards TM7, thereby strengthening interactions 
between TM6 and TM7 and between TM7 and TM2 (Figure 2C, 
D). In contrast, interactions between TM3 and TMS, measured as 
a hydrogen bond between Thr^'^' and Ser^ *^, remained mostiy 
unmodified (Figure 2C). 

The predicted binding modes for the endogenous agonist 
dopamine and the two partial agonists aripiprazole and FAUC3S0 
differ significantiy from that of the antagonist eticlopride in the 
crystal structure of the closely related dopamine D3 receptor (D3R) 
[35] (Figure S8). The fuU agonist dopamine is located deeper 
within the binding pocket of D2R than eticlopride (Figure S8A). In 
this conformation, the ammonium nitrogen of dopamine is 
stabilized by an ionic interaction to Asp^'^^ and a cation-jt 
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Figure 1. Ligand-specific conformations of the extracellular surface above the binding pocket. Representative snapshots of the 
dopamine- (left), the aripiprazole- (middle) and the FAUC350-complexes (right) are shown (average structures taken from 950-975 ns, 750-775 ns 
and 350-375 ns, respectively). (A) Extracellular view from the top into the binding pocket of the simulation systems. The ligands dopamine, 
aripiprazole and FAUC350 are highlighted in an orange, dark-grey and blue balls and sticks mode, respectively. The backbone of D2R is shown as 
ribbons, with important amino acids stabilizing the ligands indicated as sticks. (B) Side view into the binding pocket of the simulation systems, with a 
longitudinal section through the surface of D2R. Dopamine, aripiprazole and FAUC350 are represented as orange, dark-grey and blue balls and sticks, 
respectively. Compared to dopamine, both aripiprazole and FAUC350 open up the binding pocket towards the extracellular surface. 
doi:10.1371/journal.pone.0100069.g001 
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Figure 2. Investigation of dihedral angles and hydrogen-bond networks within the ligand binding pocket of the simulation 
systems. (A) Ligand-specific regulation of the dihedral angle of residue His393^"^^ (atoms: C-Ca-Cp-Cy), depicted as red, grey and blue lines for the 
dopamine-, the aripiprazole- and the FAUC350-complexes, respectively. (B) The dihedral angle of residue Serl94^'*^ (atoms: C-Cot-Cp-Cy) within the 
dopamine-simulation is shown as red lines. Unlike a constant value observed within the dopamine system, both aripiprazole (left insert) and FAUC350 
(right insert) cause a greater flexibility of this dihedral angle. (C) Hydrogen-bond interactions between representative residues of helices TIVI2, TiV13, 
TM5, TIVI6 and TM7. Aripiprazole (grey values) and FAUC350 (blue values) cause a ligand-specific modulation within interhelical networks in proximity 
to the binding pocket compared to dopamine (red values). (D) Representative snapshots of D2R within the dopamine-, the aripiprazole- and the 
FAUC350-complexes shown as red, grey and blue ribbons, respectively. The snapshots represent average structures taken from 950-975 ns 
(dopamine), 750-775 ns (aripiprazole) and 350-375 ns (FAUC350). The superposition of these structures visualizes ligand-specific changes within 
interhelical networks in proximity to the binding pocket. Helix movements are indicated with green arrows. 
doi:1 0.1 371/journal.pone.01 00069.g002 
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interaction to Phe^'^', thereby connecting TM3 and TM6 at 
residues that have been shown to be important for both ligand 
binding and receptor activation [24,36]. This connection between 
two helices that are relevant for activation is prevented by 
eticlopride because its positively charged nitrogen atom is 
displaced by 3.4 A compared its dopamine equivalent and its 
ethyl moiety helps shield the positive charge. Comparable but less 
pronounced displacements of the positively charged nitrogen 
atoms were observed in the partial agonist complexes (2.2 A and 
2.9 A for aripiprazole and FAUC350, respectively) (Figure S8B 
and C). In these cases, it is still possible to connect helices TM3 
and TM6 structurally, even though the conformation of Phe^'^' 
differs from that within the dopamine complex (and His^'^^ helps 
stabilize the charged nitrogen atoms by cation- Jt interactions). 
These observations thus suggest a more general role for the 
positively charged nitrogen atom of (partial) agonists in stabilizing 
the active-state of D2R in contrast to its role in the binding mode 
of antagonists such as eticlopride. Further studies of N-substituted 
agonists are needed in order to explore the possible conforma- 
tional influences of the structural connections around the 
protonated nitrogen on receptor activation. 

A Hydrophobic Network Connects the Ligand- and the G 
Protein-binding Poclcets and is Regulated Differently by 
the Ligands 

Ligand-binding to the extracellular part of the receptor is 
connected to conformational changes on its intracellular side [37], 
which points to the existence of an allosteric communication path 
that transforms rather small changes within the orthosteric binding 
pocket into pronounced intracellular rearrangements [38-41]. 
Earher studies have identified hydrophobic residues at the core of 
TM3, TM5, TM6 and TM7 to be involved in this signal 
propagation; key roles in receptor activation were attributed to the 
so called 'transmission switch', consisting of He"* *", Pro^'^*^ and 
Phe^**[39], and the 'rotamer toggle switch', centered around 
Trp*^ *" [42,43]. C ()nsist(mt with these studies, our MD simulations 
depicted an allosteric communication network that links the 
Ugand-binding pocket to the G protein-coupling domain 
(Figure 3A, B). Starting from distinct dihedral angles of His^'^^ 
within the binding pocket, we observed ligand-specific conforma- 
tional changes of individual residues of this network (Figure S9), 
including the aromatic amino acids Phe ' , Trp ' and Phe' ' , 
which are known to be crucial for receptor a(ti\ation 
[39,40,44,45]. The lower end of this network is formed by the 
highly conserved residues Tyr^'^^ and Tyr''^^ (Figure 3C), which 
were suggested to stabilize the active-state of the receptor via a 
water-mediated hydrogen-bond [46,47] . In analogy, it was found 
that mutation of Tyr ' in rhodopsin is involved in allosteric 
coupling to EL2[48] and in a reduction in the capacity to activate 
transducing [49]. Whereas the hydrogen bond between Tyr^'^^ 
and Tyr''^^ remained stable throughout the dopamine-simulation, 
aripiprazole and FAUC350 caused a larger fluctuation in the 
distance between these residues (Figure 3D). 

Full and Partial Agonists Influence the Conformation of G 

Protein-binding Epitopes of D2R Differently 

Finally, ligand-specific conformational changes involve domains 
of receptor-G protein coupling. The crystal structure of a 
representative ternary signaling complex [21] provides us with a 
precise molecular understanding of the interactions between an 
activated receptor and its G protein. Numerous experimental 
studies indicated that some receptor-domains constitute critical 
determinants for the activation of G proteins, including intracel- 



lular loop 2 (IIj2) [50], intracellular loop 3 (IL3, connected to the 
intracellular ends of TM5 and TM6) [51-54], the DRY-motif 
(located at the intracellular end of TM3) [55,56] and the proximal 
part of helix 8 (H8) [57]. Consequently, we analyzed the impact of 
our ligands on the conformation of these domains and found that 
the partial agonists aripiprazole and FAUC350 induce similar 
receptor conformations. However, these structures differ signifi- 
cantly from the receptor-dopamine complex (Figures 4 and 5A). 
One of these differences refers to the conformation of Metl40''"^, 
whose side chain formed contacts to the G protein in the 
dopamine simulation, but was directed away in the aripiprazole 
and FAUC350 complexes (Figures 4A-D and SIO). Moreover, a 
computational alanine scanning analysis of this residue revealed an 
impaired stabilization of the receptor-G protein interface within 
the dopamine-complex, whereas the Ml 40 A mutation within the 
aripiprazole- and the FAUC350-complexes exhibited weaker 
effects on receptor-G protein coupling (Figure SIO), indicating a 
less important role for Metl40'^^ in stabilizing the receptor-G 
protein interface within the latter two simulation systems. These 
observations are consistent with experimental studies showing that 
mutation to alanine at the corresponding position of f)2AR and the 
muscarinic receptors Ml and M3, resulting in a loss of interaction, 
was connected to a reduced capacity to active G proteins [50] . IL2 
is coupled structurally to the highly conserved DRY-motif at the 
intracellular end of TM3 (Figure 4A-C), which is known to be 
involved in G protein coupling via Arg'^ '^" [55,56]. In the case of 
the dopamine and aripiprazole complexes, although with a subdy 
rearranged architecture in the aripiprazole simulation due to a 
conformational change of IL2 around MetHO''^^, Arg^'*'' was 
stabilized by Asp'* and Tyrl42'^^ and formed an ionic 
interaction to Asp350 of the C-terminal part of Gcx; (Figure 4A- 
C, Figure Sll). This stabilizing triad was mainly broken in the 
FAUC350-simulation, resulting in a less stable salt bridge (Figure 
SI 1). In addition to the intracellular part of TM3, the C-terminus 
of Got; is surrounded by the intracellular ends of TM5 and TM6 
(constituting the beginning and the end of 1L3, respectively) and 
the junction of TM7 and H8. Contacts of the C-terminus of Got; to 
the intracellular domain of TM5 were maintained throughout all 
three simulations (but with conformational changes for the open 
end of the N-terminal part of IL3) (Figures 4D and SI 2). However, 
the junction of TM7 and H8 moved away from the C-terminus of 
Ga; in the aripiprazole and the FAUC350 complexes (Figure 4F). 
Moreover, cation-7t (Lys ' /Phe354) and hydrophobic interac- 
tions (Met^"'^) between residues of TM6 and the C-terminal part 
of Got; were reduced in the aripiprazole- and the FAUC350- 
binding models (Figure 4A— C and SI 3). 

Taken together, the conformational changes described in the G 
protein-binding domains of D2R directly influenced the confor- 
mation of the interacting epitopes of Gaj, including the C- 
terminus, helices a4 and a5 and the loops ocN/pi and P2/p3 
(Figure 4). As it was recently shown that receptor-catalyzed 
nucleotide exchange is transmitted via dynamic changes within the 
linker regions connecting the areas of receptor-G protein and 
nucleotide-G protein-coupling [3,5], we hypothesize that a 
complete G protein activation requires specific intracellular 
receptor conformations, which can only be stabilized by a fuU 
agonist like dopamine (Figure 5A). Distinct partial agonist-induced 
differences in the way intracellular epitopes are shaped may lead to 
an impaired receptor-G protein coupling and thus mcidulate the 
extent of the functional response (Figure 5B). It is thcrc-fore 
tempting to assume that impaired receptor-G protein coupling due 
to distinct ligand-specific conformations is a major determinant of 
partial agonist efficacy. 
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Figure 3. Hydrophobic network between the extracellular and the intracellular surface of D2R. The snapshots shown represent average 
structures taken from 950-975 ns (dopamine), 750-775 ns (aripiprazole) and 350-375 ns (FAUC350). (A, B) The backbone of D2R is shown as ribbons, 
important amino acids comprising the hydrophobic network are visualized as sticks. The ligands within their binding pockets are highlighted as 
dotted spheres. Superposition of representative snapshots taken from the dopamine- (red), the aripiprazole- (grey) and the FAUC350-simulations 
(blue) indicate ligand-specific conformations of residues within this hydrophobic network. (C) A water-mediated hydrogen bond between residues 
Tyr^^** and Tyr"^ (represented as sticks) of the crystal structure of P2AR bound to the ligand BI167107 and an intracellular nanobody (PDB-ID: 4LDE) 
is shown in green, with residue Arg^"^° forming the upper end of the G protein binding pocket. Additionally, representative snapshots of the 
dopamine- (red), the aripiprazole- (grey) and the FAUC350-simulations (blue) are superimposed with the crystal structure. (D) The distance between 
the hydroxyl groups of residues Tyr"^ and Tyr"^ within the dopamine-, the aripiprazole- and the FAUC350-complexes are shown as red, grey and 
blue lines, respectively. 
doi:1 0.1 371 /journal.pone.01 00069.g003 




Figure 4. Crucial amino-acid interactions and conformational changes within G protein coupling domains. The backbone of D2R is 
shown as ribbons, whereas important amino acids are highlighted as sticks. The snapshots provided represent average structures taken from 950- 
975 ns (dopamine), 750-775 ns (aripiprazole) and 350-375 ns (FAUC350). (A-C) Crucial interactions between residues from the C-terminal part of Got 
(Asp350, Phe354) and residues from the DRY-motif of TM3 (Asp^"', Arg^^°), from IL2 (IVIetMO, Tyr142) and from TIVI6 (Lys"^ Met"*") within 
representative snapshots of dopamine-, aripiprazole- and FAUC350-complexes are depicted in red, grey and blue, respectively. (D) A superposition of 
representative snapshots of dopamine-, aripiprazole- and FAUC350-complexes, represented in red, grey and blue, respectively, indicates 
conformational changes for residue MetMO of IL2 and for the N-terminal part of IL3. (E) Enlarged view of Figure 4d on the conformational changes of 
residue IVIetMO of IL2 within the simulation systems. A green arrow visualizes the movement of residue IVIetMO. (F) A superposition of representative 
snapshots of dopamine-, aripiprazole- and FAUC350-complexes, represented in red, grey and blue, respectively, highlights the increasing distance of 
the intracellular part of TM3 and the junction of TM7 and H8, measured as the distance between the Cot-atoms of Ala135^'^^ and Ile431 of TM7/H8 
(dashed box). 

doi:1 0.1 371/journal.pone.01 00069.g004 
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Figure 5. Conformational changes within G protein coupling domains and their supposed effect on nucleotide release. (A) A 

summary of certain conformations of G protein coupling epitopes of D2R observed within the simulation systems is shown. (B-D) One representative 
snapshot of the dopamine-complex, taken as an average structure from between 950-975ns, is used as a scaffold, in which to compare the effect of 
the ligands dopamine, ariplprazole and FAUC350 on the conformation of G protein coupling domains and thus on nucleotide release schematically. 
Areas of receptor-G protein coupling are shown as dark-grey and light-grey ribbons, respectively. The conformation of GDP has been taken from a 
superposition of the aforementioned snapshot with the crystal structure of ground-state Goti (PDB-ID: 1GP2). Stable contact regions of D2R are 
highlighted in green, conformational changes are indicated in orange. Colored arrows imply the supposed contribution of individual G protein 
coupling domains (contacts to cxS: red, contacts to czN/pl: yellow, contacts to a4/(56: violet) on nucleotide release (green arrows). 
doi:1 0.1 371 /journal.pone.01 00069.g005 



Conclusions 

In summary, we have used representative homology models of 
ternary receptor-G protein-complexes as structural scaffolds to 
investigate the molecular basis of partial agonism. We were able to 
capture distinct ligand-specific conformations within a homology 
model of our recendy described ternary D2R-Gai-complex, which 
help explain the graded efficacy of 1,4-DAP partial agonists such 
as ariplprazole and FAUC350 in comparison to the fuU agonist 
dopamine. 

However, ligand-induced structural changes may differ for 
other receptor-effector systems such as receptor- (3-arrestin-com- 
plexes [58] or even complexes of receptors with other G protein 
subtypes. This consideration is also relevant for both ariplprazole 
and FAUC350, which have been shown to exhibit biased signaling 
properties at D2R with respect to the activation of G protein- or P- 
arrestin-pathways, to Gai/Got„-signaling, or the stimulation of 
ERKl/2 phosphorylation [25,59,60]. Therefore, future studies 
will be required, ideally based on atomistic templates, to sample 
the conformations of a certain receptor-effector-complex entirely 
in order to explore the molecular determinants of biased agonism. 



One purpose of our work is to investigate the potential of long 
MD simulations for investigating complex biological processes, 
especially for GPCRs, for which experimental evidence is often 
sketchy. The simulations reported above are at the high end of 
what is possible today, appear inherendy reasonable and offer 
rationalizations of experimental observations. We have concen- 
trated on "hard" results (structures, persistent interactions) in the 
main text and have reported less well-founded data (e.g. 
MMPBSA results) in the Supporting Information in order to 
provide as reliable results as possible. However, the simulations are 
inherendy stochastic (because of their starting velocities) and the 
force fields only moderately well tested for simulations of this 
length. In particular, our preferred strategy of using single long 
simulations is not without alternative. 

Nevertheless, we believe that the simulations reported above are 
relevant for the real GPCR system and that they potentially 
provide new atomistic details that can now be tested experimen- 
tally. 
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Materials and Methods 

Simulation Systems 

The simulation systems contain our recently described active- 
state homology models of the dopamine D2 receptor (D2'^"""R 
and D2^''R, depending on the initial rotamer conformation of the 
side chain of residue His393'''^^ in the D2R models) in complex 
with a nucleotide-free Gofij-protein [22], which were based on the 
crystal structure of the P2-adrenergic receptor (P2AR) together 
with a heterotrimeric G^-protein [21]. In addition, the models are 
embedded in a hydrated membrane consisting of dioleoylpho- 
sphatidylcholine (DOPC) lipids and coupled to the ligands 
dopamine, aripiprazole or FAUC350 (Tables SI and S2). 

Simulation A and B (Table SI) refer to previously published 
simulations of 1 ^s each [22]. For simulation A, we performed an 
additional 500 ns simulation run as described [22]. 

Simulation systems C and D (Table SI) were prepared as 
follows: The ligands aripiprazole and FAUC350 were geometry 
optimized by means of Gaussian 09 [61] at the HF/ 6-3 1 (d,p) level 
(attributing a formal charge of + 1). AutoDock Vina [62] was used 
to subsequently dock both ligands into a membrane-inserted 
conformation of simulation system B. The ligand dopamine was 
removed before the docking procedure. We applied a search space 
of 28x26x40 A to ensure a complete coverage of the binding 
pocket. The ligands were subjected to the docking procedure using 
an exhaustiveness value of 32 and a randomly selected starting 
position. 20 conformations of each ligand were obtained and 
inspected manually. Based on the scoring function of AutoDock 
Vina and experimental data, we selected one final conformation 
for each ligand. Parameter topology and coordinate files for the 
docked complexes were build up using the tleap module of 
AMBER 10 [63] and subsequently converted into GROMACS 
input files [64,65] . We finally exchanged the coordinates of the 
ternary dopamine-D2^''R-Gaii-complexes (system B) within the 
membrane-inserted simulation systems with those of the docked 
aripiprazole- and FAUC350-D2"PR-Gan-compl exes. The final 
simulation systems contained 460 DOPC-lipids surrounding the 
proteins and 8 chlorine atoms for charge neutralization. In total, 
systems C and D consisted of 227,577 atoms (51,300 water 
molecules) and 227,571 atoms (51,298 water molecules), respec- 
tively. 

The final simulation systems were submitted to energy 
minimization (2500 steps of steepest descent minimization), 
equilibration (10 ns) and production molecular-dynamics simula- 
tion runs (800 ns and 500 ns for system C and D, respectively) 
using the GROMACS simulation package [64] as described 
earlier [26]. For all simulations, the general AMBER force field 
(GAFF) [66] was used for the ligands and the DOPC molecules 
and the force field f[99SB [67] for the protein residues. The GAFF 
force field for the lipids has been validated extensively by the 
original authors [68] and in our earlier work [22,26]. The SPC/E 
water model [69] was appUed. Parameters for the ligands were 
assigned using antechamber [63] and charges were calculated 
using Gaussian 09 [61] at the HF/6-31(d,p) level and the RESP 
procedure according to the literature [70] . A formal charge of + 1 
was defined for the ligands. Throughout the productive simula- 
tions, a force of 1.0 kcal mol~' was applied to the N-terminal 
part of the G-protein's ocN-heUx as described previously [22]. 

Data Analysis 

We removed water and DOPC molecules for data analysis. The 
analysis of the trajectories was performed with the PTRAJ module 
of AMBERIO [63]. Calculation of the binding free energies was 



accompUshed using MMPBSA. Py [71]. Figures were prepared 
using PyMOL [72] and Chimera [73]. 

Supporting information 

Figure SI Analysis of the dopamine simulations A and 

B. (A-B) Representative conformations of the binding pocket of 
D2R within the simulation systems A and B are shown in green 
and red, respectively. Residues stabilizing dopamine (shown as 
sticks) in its binding pocket are represented as sticks, the backbone 
of D2R is shown as ribbons. Whereas both hydroxyl groups of 
dopamine's catechol moiety participate in stabilizing a hydrogen 
bond network comprised of residues Ser^ *^, Ser'^'^, Ser' **' and 
His^"'''' in simulation B, dopamine is forming only one stable 
hydrogc-n bond to residue Ser^ *^ in simulation A. (C) A hydrogen- 
bond analysis between dopamine and residues occupying the 
binding pocket of D2R is provided. (D) The dihedral angle of 
residue His393''" (atoms: C-Ca-CfS-Cy) for simulation A and B is 
depicted as green and red lines, respe[:tively. (E— F) Representative 
conformations of the intracellular part of D2R within the 
simulation systems A and B are shown in green and red, 
respectively. Important amino acids are visualized as sticks. A 
(water-mediated) hydrogen bond between residues Tyr^'^® and 
Tyr''^^ of D2R and a salt bridge between residue Arg^' " of D2R 
and Asp350 of Ga was only observed within simulation B, but not 
within simulation A. (G) The distances between the hydroxyl 
groups of the tyrosines Tyr '" ^** and Tyr'" '"' of D2R are depicted as 
green and red Knes, respectively. (H-I) Free energy of binding 
calculations have been performed for dopamine-D2R (H) and 
D2R-Gai (I) using the GBSA-Method. The values are shown as 
green and red lines for simulation A and B, respectively, and 
indicate, in both cases, more stable interactions within simulation 
B. 

(TIFF) 

Figure S2 RMS-deviations of the simulation systems. 

RMS-deviations for individual components of the simulation 
systems are shown. The ligands and the receptors are fitted on the 
Ca-atoms of the receptors, whereas the G proteins are fitted on the 
Ca-atoms of the G-proteins. (A) RMSD-values for the ligand 
dopamine, D2R and Got; are given in yellow, dark-red and 
salmon, respectively. (B) RMSD-values for the ligand aripiprazole, 
D2R and Ga; are given in black, dark-grey and light-grey, 
respectively. (C) RMSD-values for the ligand FAUC350, D2R and 
Goc; are given in turquoise, dark-blue and light-blue, respectively. 
The values for FAUC350 indicate the existence of two 
interconvertible ligand conformations. 
(TIFF) 

Figure S3 Atomic fluctuations within the simulation 
systems. Atomic fluctuations for the Ca-atoms of the dopamine- 

(A), the aripiprazole- (B) and the FAUC350-complex (C) are 
shown in red, grey and blue, respectively. The thick lines for 
receptors and G proteins refer to a fit on Ca-atoms of receptors 
and G proteins, respectively, whereas the thin lines represent the 
fluctuations of the G proteins fitted on the receptor moieties. 
(TIFF) 

Figure S4 Distances between receptors and the G- 
termini of the G proteins. Distances between the centers of 

mass of D2R and the C-terminus of Ga; for the dopamine- (A), the 
aripiprazole- (B) and the FAUC350-complex (C) are shown in red, 
grey and blue, respectively. 
(TIFF) 

Figure S5 The binding pocket of the simulation sys- 
tems. Side view into the binding pocket of the simulation systems. 
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The backbone of D2R is shown as ribbons, important amino-acids 
stabilizing the conformation of the ligands are represented as 
sticks. (A) A representative snapshot of the conformation of 
dopamine (orange balls and sticks) within the binding pocket is 
shown. (B) A representative snapshot of the conformation of 
aripiprazole (dark-grey balls and sticks) within the binding pocket 
is visualized. (C, D) Representative snapshots of the conformation 
of FAUC350 within the binding pocket are highlighted, taken 
from within the last 20 ns of the simulation time (light-blue balls 
and sticks, C) and at 400 ns (blue balls and sticks, D). 
(TIFF) 

Figure S6 Residues within the binding pocket of D2R 
interacting with the ligands. A detailed contact analysis of 
residues within the binding pocket of the simulation systems 
interacting with the ligands is provided. An amino acid is 
considered as forming a contact to a ligand when at least one 
atom of the amino acid approaches at least one atom of the ligand 
closer than 3.5 A. The contacts are investigated throughout the 
simulated time scales. 
(TIFF) 

Figure S7 Distance between EL2 and the extracellular 
end of TM7. The distances between the side chains of He 183 of 
EL2 and Tyr408' ' of TM7 for the dopamine-, the aripiprazole- 
and the FAUC350-complex are shown in red, grey and blue, 
respectively. 
(TIFF) 

Figure S8 Comparison of the predicted binding modes 
of our agonists at D2R with the conformation of the 
antagonist eticlopride at D3R. Side view into representative 
snapshots of the binding pockets of D2R and the crystal structure 
of D3R. The snapshots represent average structures taken from 
950-975 ns (dopamine), 750-775 ns (aripiprazole) and 350- 
375 ns (FAUC350). The backbone of the receptors is shown as 
ribbons, the ligands and important amino acids (Asp^'^^, Phe^'^' 
and His^"^^) stabilizing their conformation are represented as 
sticks. The positively charged nitrogen atoms of the ligands are 
highlighted as blue balls, wh(;rcas th(^ distances between these 
nitrogen atoms of the ligands arc gi\ i'n in light pink. The figure 
shows an overlay of eticlopride (green) at D3R with dopamine 
(orange/red. A), aripiprazole (grey sticks, B) and FAUG350 (blue, 
C) at D2R. 
(TIFF) 

Figure S9 Ligand-specific dihedral angles of represen- 
tative residues comprising the hydrophobic network 
between the ligand and the G protein binding pockets. 

(A-H) Dihedral angles (atoms: C-Ca-CP-Cy) of important 
residues from the core of the hydrophobic network (lle'^*", 
Tyr^ «, Phe^ **, Tip^ «, Phe^"'^', Phe^^ ''^ Tyr' '^^ Tyr^ '^'^ which 
connect the ligand and the G protein binding pockets of D2R, are 
shown as dark-red, dark-grey and dark-blue lines representing the 
dopamine-, the aripiprazole- and the FAUC350-complexes, 
respectively. In addition, the dihedral angle of residue Trp*' *'' 
between atoms Ca-CP-Cy-CSj (D) is provided as light-red, light- 
grey and light-blue lines for the dopamine-, the aripiprazole- and 
the FAUC350-complexes, respectively. 
(TIFF) 

Figure SIO Investigations on residue Metl40 of IL2. (A) 

A computational alanine scanning analysis for residue Met 140 of 
IL2 is provided for the dopamine- (left), the aripiprazole- (middle) 
and the FAUC350-complexes (right). Whereas the M140A 



mutation within the dopamine-complex was connected to an 
impaired stabilization of the receptor-G protein interface, we 
observed weaker effects of this mutation within the aripiprazole- 

and the FAUC350-complexes indicating a less important role of 
Metl40""^ within the latter two simulation systems. (B) The 
distance between the side chain of Met 140 of IL2 and the Ca- 
atom of Ile343 of a5 for the dopamine- (left), the aripiprazole- 
(middle) and the FAUC350-complexes (right) is shown. An 
increasing distance between these residues indicates a conforma- 
tion of Met 140 exhibiting reduced contacts towards the G protein. 
(TIFF) 

Figure SIl Distances between individual residues of 
TM3, IL2 and the a5 helix of the G protein. The distances 
between individual residues of TM3 (Asp^ *^ Arg^-*°), IL2 
(Tyrl42) and the a5 helix of the G protein (Asp350) are shown 
as red, grey and blue lines for the dopamine-, the aripiprazole- and 
the FAUC350-complexes, respectively. The distance between 
residues Arg' ' ' and Asp350 of the G protein comprising a salt 
bridge (A), between residues Arg " ' and Asp " of TM3 (B), 
between residues Arg^'^" of TM3 and Tyrl42 of IL2 (C) and 
between residues Asp^ *^ of TM3 and Tyrl42 of IL2 (D) is 
provided. 
(TIFF) 

Figure S12 RMS-deviations and contact analysis of the 
TM5-IL3 region. (A) RMS-deviations for TM5 of D2R within 
the simulation .systems are shown as red, grey and blue lines for the 
dopamine-, the aripiprazole- and the FAUC350-complexes, 
respectively. The values attribute a low conformational flexibiUty 
to TM5. (B) RMS-deviations for the proximal part of IL3 of D2R 
within the simulation systems are shown as red, grey and blue lines 
for the dopamine-, the aripiprazole- and the FAUC35()-complex- 
es, respectively. The values indicate a high conformational 
flexibility for IL3. (G, D) A detailed contact analysis between 
residues of TM5 and IL3 of D2R interacting with residues of the 
G protein is provided. An amino acid is considered as forming a 
contact when at least one atom of one amino acid approaches at 
least one atom of a second amino acid closer than 3.5 A. The 
contacts are investigated throughout the simulated time scales. 
(TIFF) 

Figure S13 Investigation of TM6 residues Lys367*"'^' 
and Met374''-'*. (A) The distances between TM6 residue 
Lys367''-^^ of D2R and the C-terminal residue Phe354 of Ga 

are shown. (B) The distances between D2R residues Argl32'' '^'' 
and Met374' ' ^ are shown. Values are represented in red, grey 
and blue for the dopamine-, the aripiprazole- and the FAUC350- 
complex, respectively. 
(TIFF) 

Table SI Overview of the simulation systems and their 
simulated time scales. 

(DOCX) 

Table S2 Chemical structures of the ligands investigat- 
ed. 

(DOCX) 
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